Columbus' egg: Oral drugs are absorbed in finite time.

The infinite time of oral drug absorption was conceived from the first day of the birth of pharmacokinetics when H. Dost introduced the term pharmacokinetics in his book published in 1953. He adopted the function developed by H. Bateman back in 1908 for the decay of the nuclei isotopes to describe oral drug absorption as a first-order process. We unveiled this false hypothesis relying on common wisdom i.e. drugs are absorbed in finite time. This false assumption had dramatic effects on the evolution of oral pharmacokinetics but most importantly on the bioavailability and bioequivalence concepts and metrics. This work focuses on the finite absorption time (FAT) concept, the relevant Physiologically Based Finite Time (PBFTPK) models developed and their applications in oral pharmacokinetics, bioavailability and bioequivalence. The crux of the matter is that drug absorption from the gastrointestinal tract takes place under sink conditions because of the high blood flow rate in the vena cava. The termination of oral, pulmonary and intranasal drug absorption at a specific time point, calls for regulatory changes in bioavailability and bioequivalence studies in terms of the study design and metrics used for the bioequivalence assessment.

Enhancement of Docetaxel Absorption Using Ritonavir in an Oral Milk-Based Formulation.

OBJECTIVE: The current study aimed to develop a novel milk-based formulation of docetaxel, a sparingly soluble antineoplastic agent, administered so far exclusively by the intravenous route and evaluate its oral bioavailability. METHODS: Pre-formulation studies included the determination of docetaxel solubility in water-alcohol mixtures as well as short-term content uniformity experiments of the final formulation. The pharmacokinetic (PK) performance of the developed milk-based formulations was further evaluated in vivo in mice using ritonavir, a potent P-glycoprotein inhibitor, as an absorption enhancer of docetaxel and the marketed intravenous docetaxel formulation, Taxotere®, as a control. RESULTS: In vivo PK results in mice showed that all the administered oral docetaxel formulations had limited absorption in the absence of ritonavir. On the contrary, ritonavir co-administration given as pre-treatment significantly enhanced oral bioavailability of both the marketed and milk-based docetaxel formulations; an even more marked increase in drug exposure was observed when ritonavir was incorporated within the docetaxel milk-based formulation. The fixed-dose combination also showed a more prolonged absorption of the drug compared to separate administrations. CONCLUSIONS: The current study provides insights for the discovery of a novel milk-based formulation that could potentially serve as an alternative, non-toxic and patient-friendly carrier for an acceptable docetaxel oral chemotherapy.

On the pharmacokinetics of two inhaled budesonide/formoterol combinations in asthma patients using modeling approaches.

Dry powder inhalers containing the budesonide/formoterol combination have currently a well-established position among other inhaled products. Even though their efficacy mainly depends on the local concentrations of the drug they deliver within the lungs, their safety profile is directly related to their total systemic exposure. The aim of the present investigation was to explore the absorption and disposition kinetics of the budesonide/formoterol combination delivered via two different dry powder inhalers in asthma patients. Plasma concentration-time data were obtained from a single-dose, crossover bioequivalence study in asthma patients. Non-compartmental and population compartmental approaches were applied to the available datasets. The non-compartmental analysis allowed for an initial characterization of the primary pharmacokinetic (PK) parameters of the two inhaled drugs and subsequently the bioequivalence assessment of the two different dry powder inhalers. The population pharmacokinetic analysis further explored the complex absorption and disposition characteristics of the two drugs. In case of inhaled FOR, a five-compartment PK model including an enterohepatic re-circulation process was developed. For inhaled BUD, the incorporation of two parallel first-order absorption rate constants (fast and slow) for lung absorption in a two-compartment PK model emphasized the importance of pulmonary anatomical features and underlying physiological processes during model development. The role of potential covariates on the variability of the PK parameters was also investigated.

Population pharmacokinetics of fluticasone propionate/salmeterol using two different dry powder inhalers.

The combination of fluticasone propionate (FLP) and salmeterol (SAL) is often used in clinical practice for the treatment of pulmonary disorders. The purpose of this study was to explore the pharmacokinetics (PK) of inhaled FLP and SAL, after concomitant administration, in healthy male and female subjects using two dry powder inhalers. Plasma concentration (C)-time (t) data were obtained from a single dose, two-sequence, two-period, crossover (2×2) bioequivalence (BE) study. Activated charcoal was co-administered in order to prohibit absorption from the gastrointestinal tract. A number of 60 subjects were recruited, while 57 of them completed the study and were included in the PK analysis. Initially, PK parameters of FLP and SAL were estimated using the classic non-compartmental methods. Subsequently, BE assessment was applied to the estimated PK parameters of the two dry powder inhalers. Special focus was placed on the population PK analysis of the C-t data, which were pooled together. 'Treatment' (i.e., test or reference) and 'period' of the BE study were considered as covariates. A variety of structural and residual error models were tested to find the one which best described the plasma C-t data of FLP and SAL. Demographic data were also evaluated for their impact on the PK parameters. Several goodness-of-fit criteria were utilized. The non-compartmental PK estimates of this study were in agreement with previously reported values. The population PK analysis showed that FLP data were described by a two-compartment model with first-order absorption and elimination kinetics. Body weight was found to affect significantly absorption rate constant, inter-compartmental clearance, and volume of distribution of the peripheral compartment. As body weight increases, the values of these PK parameters also rise. For SAL, the best results were obtained when a two-compartment disposition model was used assuming very rapid absorption kinetics (like intravenous bolus) and first-order elimination kinetics. Gender was found to be a significant covariate on clearance, with men exhibiting higher clearance than women.

Stability and physicochemical characterization of novel milk-based oral formulations.

PURPOSE: The purpose of this work was to assess the colloidal stability of novel milk-based formulations. METHODS: Milk-based formulations were prepared in situ by adding into milk alkaline- or ethanolic-drug solutions containing an array of drugs namely; ketoprofen, tolfenamic acid, meloxicam, tenoxicam and nimesulide, mefenamic acid, cyclosporine A, danazol and clopidogrel besylate. The produced formulations were characterized by means of dynamic lightscattering, ζ-potential studies, atomic force microscopy, fluorescence spectroscopy, Raman spectroscopy complemented with ab initio calculations and stability studies. RESULTS: The presence of the drugs did not induce significant changes in most cases to the particle size and ζ-potential values of the emulsions pointing to the colloidal stability of these formulations. Raman spectroscopy studies revealed interactions of the drugs and the milk at the intermolecular level. Complementary analysis with ab initio calculations confirmed the experimental observations obtained by Raman spectroscopy. Finally the produced drug containing alkaline/ethanolic solutions exhibited stability over a period of up to 12 months. CONCLUSIONS: The current data demonstrate that milk is a promising drug carrier.

Computational-Regulatory Developments in the Prediction of Oral Drug Absorption.

Early prediction of human intestinal absorption is important in selection of potential orally administered drugs. Various computational models for prediction of the fraction of dose absorbed, Fa, have been developed. In 1989, a sigmoidal relationship between Fa and drug absorption potential was shown. Since then various physicochemical descriptors of molecules (lipophilicity, polar surface area, hydrogen bond descriptors) have been found to correlate with human intestinal absorption and various attempts in estimating Fa have been reported. Most studies rely on the presupposition that Fa is mainly dependent on drug's solubility, which drives the dissolution rate in the gastrointestinal (GI) fluids, and the rate of passive drug transport across the intestinal membrane. In the same vein, the biopharmaceutics classification system (BCS) and the relevant FDA guideline classify drugs in four categories according to their aqueous solubility and permeability. However, the biopharmaceutics drug disposition classification system (BDDCS) revealed the poor predictability of permeability estimates for Fa and the major role of transporters for GI uptake of drugs. The role of solubility in the reaction limited model of dissolution and the ubiquitous presence of supersaturated solubility-dissolution phenomena in the GI lumen, call for a more physiologically relevant consideration of GI absorption.

Development of a reaction-limited model of dissolution: application to official dissolution tests experiments.

A reaction-limited model for drug dissolution is developed assuming that the reaction at the solid-liquid interface is controlling the rate of dissolution. The dissolution process is considered as a bidirectional chemical reaction of the undissolved drug species with the free solvent molecules, yielding the dissolved species of drug complex with solvent. This reaction was considered in either sink conditions, where it corresponds to the unidirectional case and the entire amount of the drug is dissolved, or reaching chemical equilibrium, which corresponds to saturation of the solution. The model equation was fitted successfully to dissolution data sets of naproxen and nitrofurantoin formulations measured in the paddle and basket apparatuses, respectively, under various experimental conditions. For comparative purposes these data were also analyzed using three functions based on the diffusion layer model. All functions failed to reveal the governing role of saturation solubility in the dissolution process associated with the diffusion layer model when the conditions for the valid estimation of saturation solubility, established theoretically in this study, were met by the experimental set up employed. Overall, the model developed provides an interesting alternative to the classic approaches of drug dissolution modeling, quantifying the case of reaction-limited dissolution of drugs.

Modelling and simulation in drug absorption processes.

Drug absorption is a complex process dependent upon drug properties such as solubility and permeability, formulation factors, and physiological variables including regional permeability differences, pH, luminal and mucosal enzymes, and intestinal motility, among others. Despite this complexity, various qualitative and quantitative approaches have been proposed for the estimation of oral drug absorption. These approaches are reviewed in this article with particular emphasis on drug dissolution modelling, dynamic systems for oral absorption and absorption models based on structure. The regulatory aspects of oral drug absorption and in particular the biopharmaceutic classification of drugs are also discussed. Models for drug dissolution and release describe adequately the in vitro data, and models for oral drug absorption provide reasonable results. The development of in vitro-in vivo correlations based on the official compendia specifications are facilitated using commercial computer packages.

In-vitro study on the competitive binding of diflunisal and uraemic toxins to serum albumin and human plasma using a potentiometric ion-probe technique.

The competitive binding of diflunisal and three well-known uraemic toxins (3-indoxyl sulfate, indole-3-acetic acid and hippuric acid) to bovine serum albumin (BSA), human serum albumin (HSA) and human plasma was studied by direct potentiometry. The method used the potentiometric drug ion-probe technique with a home-made ion sensor (electrode) selective to the drug anion. The site-oriented Scatchard model was used to describe the binding of diflunisal to BSA, HSA and human plasma, while the general competitive binding model was used to calculate the binding parameters of the three uraemic toxins to BSA. Diflunisal binding parameters, number of binding sites, n(i) and association constants for each class of binding site, K(i), were calculated in the absence and presence of uraemic toxins. Although diflunisal exhibits high binding affinity for site I of HSA and the three uraemic toxins bind primarily to site II, strong interaction was observed between the drug and the three toxins, which were found to affect the binding of diflunisal on its primary class of binding sites on both BSA and HSA molecules and on human plasma. These results are strong evidence that the decreased binding of diflunisal that occurs in uraemic plasma may not be solely attributed to the lower albumin concentration observed in many patients with renal failure. The uraemic toxins that accumulate in uraemic plasma may displace the drug from its specific binding sites on plasma proteins, resulting in increased free drug plasma concentration in uraemic patients.

The power law can describe the 'entire' drug release curve from HPMC-based matrix tablets: a hypothesis.

The purposes of this study were to (i) re-examine the relevance of Higuchi equation and the power law using both simulated and experimental release data in conjunction with the linearized, in terms of t(1/2), percent of drug release plots, (ii) demonstrate that the power law describes the entire drug release profile of several experimental data, and (iii) point out a physically based hypothesis for the successful use of power law in describing the entire drug release profile. Simulated data generated from the equation of power law were further analyzed using linear regression analysis in accord with the Higuchi equation. The analysis revealed that data generated from the equation of power law can be misinterpreted as obeying the Higuchi equation. The use of power law in describing the entire drug release curve from HPMC-based matrix tablets is validated by direct fit of power law equation to published data of other authors. A hypothesis based on the nonclassical diffusion of the solutes in the HPMC matrices is used to interpret the successful use of the power law in describing the entire release profile.

Fractal volume of drug distribution: it scales proportionally to body mass.

PURPOSE: To develop the physiologically sound concept of fractal volume of drug distribution, vf, and evaluate its utility and applicability in interspecies pharmacokinetic scaling. METHODS: Estimates for vf of various drugs in different species were obtained from the relationship: vf = (v - Vpl)(Vap - Vpl)/V + Vpl where v is the total volume of the species (equivalent to its total mass assuming a uniform density Ig/mL), Vpl is the plasma volume of the species and Vap is the conventional volume of drug distribution. This equation was also used to calculate the fractal analogs of various volume terms of drug distribution (the volume of central compartment, Vc, the steady state volume of distribution, Vss, and the volume of distribution following pseudodistribution equilibrium, Vz). The calculated fractal volumes of drug distribution were correlated with body mass of different mammalian species and allometric exponents and coefficients were determined. RESULTS: The calculated values of vf for selected drugs in humans provided meaningful and physiologically sound estimates for the distribution of drugs in the human body. For all fractal volume terms utilized, the allometric exponents were found to be either one or close to unity. The estimates of the allometric coefficients were found to be in the interval (0,1). These decimal values correspond to a fixed fraction of the fractal volume term relative to body mass in each one of the species. CONCLUSIONS: Fractal volumes of drug distribution scale proportionally to mass. This confirms the theoretically expected relationship between volume and mass in mammalian species.

Nonlinear dynamics and chaos theory: concepts and applications relevant to pharmacodynamics.

The theory of nonlinear dynamical systems (chaos theory), which deals with deterministic systems that exhibit a complicated, apparently random-looking behavior, has formed an interdisciplinary area of research and has affected almost every field of science in the last 20 years. Life sciences are one of the most applicable areas for the ideas of chaos because of the complexity of biological systems. It is widely appreciated that chaotic behavior dominates physiological systems. This is suggested by experimental studies and has also been encouraged by very successful modeling. Pharmacodynamics are very tightly associated with complex physiological processes, and the implications of this relation demand that the new approach of nonlinear dynamics should be adopted in greater extent in pharmacodynamic studies. This is necessary not only for the sake of more detailed study, but mainly because nonlinear dynamics suggest a whole new rationale, fundamentally different from the classic approach. In this work the basic principles of dynamical systems are presented and applications of nonlinear dynamics in topics relevant to drug research and especially to pharmacodynamics are reviewed. Special attention is focused on three major fields of physiological systems with great importance in pharmacotherapy, namely cardiovascular, central nervous, and endocrine systems, where tools and concepts from nonlinear dynamics have been applied.

A stochastic model describes the heterogeneous pharmacokinetics of cyclosporin.

The pharmacokinetics of cyclosporin (CsA) are unusual because of several heterogeneous features which include the presence of more than one conformer, considerable accumulation in erythrocytes and lipoproteins, extensive plasma protein binding, distribution into deep tissues, biliary secretion and hepatic clearance involving a large number of metabolites. In this study, a stochastic compartmental model was developed to describe the heterogeneous elimination kinetics of CsA. This new approach relies on a probabilistic transfer model with a gamma distributed probability intensity coefficient for drug elimination. For comparative purposes both the stochastic model and compartmental deterministic models were fitted to real post infusion data from patients receiving CsA as a 2-hr intravenous infusion. The criteria for selecting the best model showed that the stochastic model, although simpler than the compartmental deterministic models, is more flexible and gives a better fit to the kinetic data of CsA than the compartmental deterministic models. The stochastic model with a random rate intensity coefficient adequately describes the heterogeneous pharmacokinetics of CsA.

Non-linear regression analysis with errors in both variables: estimation of co-operative binding parameters.

Four different parameter estimation criteria, the geometric mean functional relationship (GMFR), the maximum likelihood (ML), the perpendicular least-squares (PLS) and the non-linear weighted least squares (WLS), were used to fit a model to the observed data when both regression variables were subject to error. Performances of these criteria were evaluated by fitting the co-operative drug-protein binding Hill model on simulated data containing errors in both variables. Six types of data were simulated with known variances. Comparison of the criteria was done by evaluating the bias, the relative standard deviation (S.D.) and the root-mean-squared error (RMSE), between estimated and true parameter values. Results show that (1) for data with correlated errors, all criteria perform poorly; in particular, the GMFR and ML criteria. For data with uncorrelated errors, all criteria perform equally well with regard to the RMSE. (2) Use of GMFR and ML lead to lower values for S.D. but higher biases compared with WLS and PLS. (3) WLS performs less well when equal dispersion is applied to the two observed variables.

Ursodeoxycholic acid modulates cyclosporin A oral absorption in liver transplant recipients.

The aim was to study the ursodeoxycholic acid (UDC) effect on the cyclosporin A (CsA) pharmacokinetics after oral administration of the microemulsion formulation Neoral (CsA-ME) in liver transplant recipients, and test the potential protective effect of this bile acid on liver and renal CsA-ME-induced toxicity. At entry into the study, 12 patients who underwent orthotopic liver transplantation received CsA-ME, for at least 6 months. They then received a cotreatment CsA-ME plus UDC (13.8 mg x kg(-1) x day(-1)) for three months. Blood concentrations of CsA were measured using a monoclonal antibody specific for the parent compound. The kinetic data were analysed by a mathematical model incorporating a time dependent rate coefficient for CsA intestinal absorption, before and after UDC treatment. Changes in serum markers of hepatic and renal injury were assessed. Individual serum bile acids were determined by chromatography. Serum levels of UDC increased from 3 to about 45% of total serum bile acids after UDC treatment. The estimated model parameters indicate that UDC administration modulates CsA intestinal absorption. In the nine non-cholestatic patients, UDC reduced the absorption rate and the bioavailability of CsA without modifying the elimination rate constant of CsA and the CsA pre-drug levels. In contrast, in the three cholestatic patients, the bioavailability tended to be higher and the absorption rate faster when CsA was combined with UDC. UDC significantly decreased elevated gamma-glutamyl transferase and creatinine serum levels and induced some clinical improvements such as disappearance of headaches in four patients. In conclusion, a 3-month UDC treatment modifies CsA intestinal absorption without affecting CsA elimination rate constant. On the other hand, UDC supplementation appears to improve CsA tolerability.

A heterogeneous tube model of intestinal drug absorption based on probabilistic concepts.

PURPOSE: To develop an approach based on computer simulations for the study of intestinal drug absorption. METHODS: The drug flow in the gastrointestinal tract was simulated with a biased random walk model in the heterogeneous tube model (Pharm. Res. 16, 87-91, 1999), while probability concepts were used to describe the dissolution and absorption processes. An amount of drug was placed into the input end of the tube and allowed to flow, dissolve and absorb along the tube. Various drugs with a diversity in dissolution and permeability characteristics were considered. The fraction of dose absorbed (Fabs) was monitored as a function of time measured in Monte Carlo steps (MCS). The absorption number An was calculated from the mean intestinal transit time and the absorption rate constant adhering to each of the drugs examined. RESULTS: A correspondence between the probability factor used to simulate drug absorption and the conventional absorption rate constant derived from the analysis of data was established. For freely soluble drugs, the estimates for Fabs derived from simulations using as an intestinal transit time 24500 MCS (equivalent to 4.5 h) were in accord with the corresponding data obtained from literature. For sparingly soluble drugs, a comparison of the normalized concentration profiles in the tube derived from the heterogeneous tube model and the classical macroscopic mass balance approach enabled the estimation of the dissolution probability factor for five drugs examined. The prediction of Fabs can be accomplished using estimates for the absorption and the dissolution probability factors. CONCLUSIONS: A fully computerized approach which describes the flow, dissolution and absorption of drug in the gastrointestinal tract in terms of probability concepts was developed. This approach can be used to predict Fabs for drugs with various solubility and permeability characteristics provided that probability factors for dissolution and absorption are available.

Modeling of supersaturated dissolution data.

A recursion equation which relies on the population growth model of dissolution is used for the analysis of supersaturated dissolution data. The concentration-time data of dissolution experiments are initially transformed to fractions of dose dissolved-generations by adopting an appropriate time interval as the time step of the recursion equation. A computer program is used to derive estimates for the maximum fraction of dose dissolved and the fraction of dose remaining in solution at steady state. Good fittings were observed when this equation was applied to phenytoin and nifedipine supersaturated dissolution data obtained from literature.